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Abstract We present a Monte-Carlo implementation 
of the Statistical Hadronization Model in e + e~ collisions. 
The physical scheme is based on the statistical hadroni- 
zation of massive clusters produced by the event gene- 
rator Herwig within the microcanonical ensemble. We 
present a preliminary comparison of several observables 
with measurements in e + e _ collisions at the Z peak. Al- 
though a fine tuning of the model parameters is not 
carried out, a general good agreement between its pre- 
dictions and data is found. 

Keywords Statistical Hadronization Model • Micro- 
canonical ensemble • Hadronization models • Hadron 
production • e + e~ annihilation • Monte-Carlo 

PACS 12.40.Ee - 05.30.-d - 13.66.Bc 



1 Introduction 

The Statistical Hadronization Model (SHM) has proved 
to be successful in reproducing essential features of the 
hadronization process, such as the multiplicities and 
the transverse momentum spectra of many hadronic 
species in high energy elementary and nuclear colli- 
sions (see e.g. [1] and references therein). Recently, it 
has been shown that this model is also able to satisfac- 
torily reproduce the rates of multi-hadronic exclusive 
channels in e + e _ annihilation at low energy [2]. Accor- 
ding to the SHM, hadronization proceeds through the 
formation of extended massive colorless objects called 
clusters or fireballs emitting hadrons according to a 
pure statistical law, i.e. all multihadronic states com- 
patible with conservation laws are equally likely. The 
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number and properties of these clusters are determined 
by the pre-hadronization dynamical process and cannot 
be predicted within the SHM. Therefore, in e + e~and 
pp collisions at high energy (y/s > 10 GeV), the SHM 
calculations have been carried out with supplementary 
assumptions, in order to obtain analytical or semi-analy- 
tical formulae to be compared directly with the data. 
Such is the case of the so-called Equivalent Global Clus- 
ter picture [U|3], in which a further assumption is in- 
voked concerning the distribution of charges and mas- 
ses among the hadronizing pre-hadronic clusters; this 
allows to make calculations in the more manageable 
canonical ensemble. 

While these supplementary assumptions have been 
very useful to point out the statistical features of ha- 
dronization, it would be desirable to perform tests of 
the SHM independently thereof, at the most funda- 
mental level of formulation. As clusters are generally 
small, calculations in a basic SHM framework involve 
the computation of averages in the microcanonical en- 
semble, i.e. with exact conservation of energy and mo- 
mentum, besides that of abelian charges (baryon num- 
ber, strangeness, electric charge, charm and beauty). 
This kind of calculation has been recently carried out 
(with the further complication of angular momentum 
conservation and more) for e + e _ annihilation at low e- 
nergy [SJ. Yet, in this case, only one cluster was in- 
volved instead of many as it is generally the case for 
high energy collisions. 

In a multi-cluster environment, the only possible 
way to cope with the problem of statistical hadroni- 
zation is a Monte-Carlo integration. In this paper a 
prototype of Monte-Carlo event generator for the SHM 
is presented and the results obtained by testing it on 
c + e~ collisions at y/s — 91.2 GeV center of mass energy, 
i.e. at the Z peak, are discussed. 
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For the present study, the developed Monte-Carlo 
hadronization code has been interfaced to the Herwig6510 
event generator [3] , which is used to simulate the e + e~ col- 
lisions except for the hadronization phase. This step is 
instead performed using the SHM hadronization mo- 
dule. This particular choice for the external generator 
is owing to the similarity between the input clusters re- 
quired by the SHM and those produced by Herwig6510 
0. In our framework these clusters are used as star- 
ting point for the statistical hadronization process. Fi- 
nally, the unstable hadron decays are performed with 
the Herwig6510 generator itself. 

The paper is organized as follows: in Sec.[2]the gene- 
ral framework of the SHM is summarized and the algo- 
rithms for the simulation of single cluster hadronization 
are presented. In Sec.[3]a comparison between the theo- 
retical formulation of the SHM and those of the Cluster 
Model implemented in other Monte Carlo generators, 
such as Herwig6510, is reported. The computational 
methods for multi-cluster hadronization in the SHM 
framework are then discussed in Sec. 0J In Sec. [5] the 
full event simulation setup, involving the external event 
generator, is described, as well as results obtained for a 
preliminary adjustment of free parameters. The impact 
on SHM predictions of variations of some Herwig6510 
free parameters involved in the cluster formation proce- 
dure is also discussed. In Sec. [6j the results obtained for 
c + c~ collisions at the Z peak with the SHM-based ha- 
dronization module are presented and compared with 
the corresponding Herwig6510 predictions and LEP ex- 
perimental data. Conclusions and future developments 
are outlined in Sec. 



2 Statistical hadronization of single clusters 

As it has been mentioned in the Introduction, the cor- 
nerstones of the SHM are: 

— in a high energy collision, as a result of the pre- 
hadronic dynamical process, a set of extended mas- 
sive colorless objects called clusters or fireballs are 
formed; 

— all multi-hadronic states confined within the cluster 
and compatible with conservation laws are equally 
likely. 

The second postulate gives rise to a definite mathe- 
matical formula [HE] involving (pseudo)projector ope- 
rators associated to conservation laws, for the proba- 
bility of observing a specific final state from a single 
cluster decay. In principle, conservation laws to be im- 
plemented include energy-momentum, angular momen- 
tum, parity and all of the internal symmetries relevant 



to strong interactions. However, like all other Monte- 
Carlo generators, we confine ourselves to abelian char- 
ges (electric charge, strangeness, baryonic number, charm 
and beauty), in addition to energy-momentum. In fact, 
in a multi-cluster environment, relevant to this work, 
dealing with non-abelian symmetry groups is overwhel- 
mingly difficult if clusters take their origin from a pure 
quantum state and are entangled with each other at the 
hadronization stage. 

In the simpler scheme of additive abelian charges 
and energy-momentum conservation, the probability of 
a multi-hadronic final state for a single cluster can be 
derived as described in [6j7]. Particularly, neglecting 
quantum statistics and quantum field effects, the mi- 
crocanonical weight ^{Nj} f° r t ne multi-hadronic chan- 
nel {Nj} made of Ni, N 2 ,-.-, N K hadrons of species 
1, . . . , K is given by [5J: 
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with P f = J2i=iPi an d Qf — Z)t=i Q») where Pi is the 
4-momentum of the i-th hadron in the channel while 
Qi = (Bi, Qi,Si, . . .) is the array collecting its abelian 
charge values (baryon number, charge, strangeness, etc.), 
and N is the total number of particles in the channel. P, 
Q and V are the cluster 4-momentum, abelian charges 
and volume respectively. Finally, it must be noted that 
a sum over particle spin states has been performed: the 
Jj quantities in the corresponding contributions to the 
total weight are the spin modulus of the single hadron 
species. The microcanonical probability p({Nj}) of the 
considered channel is given by: 
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where fi is defined as the sum over all possible channels 



(3) 



i.e. the microcanonical partition function. Finally, the 
microcanonical probability density in phase space can 
be written as: 



p(W}.M) = <4(P - Pj) ^' x 
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where {pj} = Pi,P2, ---tPn is the set of four-momenta 
of the particles in the channel {Nj}. 

The analysis of particle multiplicities has shown [T] 
that particles carrying strange valence quarks need an 
additional suppression with respect to the pure statisti- 
cal predictions. This has led to the introduction of the 
strangeness suppression parameter 7s, which, in order 
to fit its definition for inclusive multiplicities, has to 
multiply fl{Nj] of Eq. [ljas follows: 
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(5) 



where Nj is the number of hadrons of species j in the 
channel {Nj} and Sj is the total number of valence 
strange and antistrange quarks contained in the j-th 
hadron species. Extra strangeness suppression also ap- 
plies to light unflavored mesons, like the 77 meson, which 
are actually a flavorless superposition of states with a 
wave function of the general form: 



C u uu + Cddd + C s ss, 
with|C u | 2 + |C d | 2 + |CJ 2 



(6) 



1. For this kind of hadrons 
only the ss part of the wave function is supposed to be 
suppressed, with a suppression factor given by 

1 - |C S | 2 + |C S | 2 T |. (7) 

The 7s parameter and the cluster energy density p are 
the only free parameters of the model, in its microca- 
nonical formulation. The cluster proper energy density 
p plays the role of conversion factor between the cluster 
volume V, present in Eqs. l]4 
through the relationship 

M 

p = v 



and the cluster mass M 



(8) 



where Vj is the mean number of particles of the j-th 
type, K the number of included hadronic species and 
Nj the number of particles of kind j contained in the 
channel {Nj}. The mean multiplicities Vj are free pa- 
rameters which should be set with the aim of obtain- 
ing the most efficient sampling function. The abelian 
charge conservation is then introduced in the channel 
sampling procedure based on Eq.[9]in the following way: 
the whole set of hadrons is divided into eleven sampling 
groups, namely, (anti)bottomed hadrons, (anti)charmed 
hadrons, light (anti)baryons, light strange (anti)mesons, 
light charged (anti)mesons with zero strangeness and 
neutral light mesons. Moreover, the following feature 
of the multi-poissonian function is used: considering 
for simplicity only the baryon, antibaryon and meson 
hadron groups and the relative charge conservations, 
the original sampling function with an extra Kronecker 
delta, added to impose the conservation of abelian char- 
ges, can be written as 
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where 



2.1 Single cluster hadronization algorithm 

The simulation of the cluster microcanonical hadroniza- 
tion is mainly based on the algorithms described in [5] 
and here summarized. More in detail, the single cluster 
hadronization process involves two main steps, namely 
the sampling of the hadronization channel and of the 
corresponding kinematical configuration. 

Due to the large number of available decay chan- 
nels for a hadronizing cluster, an efficient sampling al- 
gorithm is needed to obtain a fast and optimized ha- 
dronization simulation: for the present case the sam- 
pling function 77 has been defined, as discussed in [5], as 
the multi-species multiplicity distribution of the SHM 
grandcanonical formulation 
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II X (N x ) = exp (-i/ x ) 



NJ. 



(11) 



is the poissonian distribution of the total number of 
baryons (x = b) or antibaryons [x = b) and where the 
corresponding mean multiplicity v x is given by the sum 
of the single baryon or antibaryon mean multiplicities. 
The functions P are the conditioned multinomial distri- 
butions of the single hadronic species, namely 



P(N?,N?,...\N x )KNJl[j^, 



(12) 



with, again, x = b for baryons and x — b for an- 
tibaryons. 

With the above decomposition of the original mul- 
tipoissonian distribution, the following sampling algo- 
rithm can be used to perform the hadronization channel 
generation: 
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1. extract randomly the number of baryons and an- 
tibaryons according to the distributions lib {Nb) and 

n-b(N-b); 

2. check whether the baryonic number is conserved. If 
not reject the sampling and go to point 1, otherwise 
generate the single baryons and antibaryons using 
the multinomial distribution of Eq. [l2| 

3. extract the single mesons using the initial multipois- 
son distributions; 

4. check the conservation of the remaining abelian char- 
ges, taking into account the already generated bary- 
ons and antibaryons. If the check is not passed start 
again with point 1, otherwise the sampled channel is 
momentarily accepted and the corresponding phase 
space availability can be verified. 

The adoption of the described sampling procedure, 
instead of the independent sampling of the single hadro- 
nic species, which would follow from Eq.[9j is motivated 
by the reduction of random number extractions which 
occur in case of rejection of the sampled channel due to 
the abelian charge conservation: the independent sam- 
pling procedure would require, for each channel sam- 
pling attempt, the extraction of K random numbers, 
where K is the number of hadron species included in 
the hadronization procedure, which would be lost in 
case of event rejection. The above algorithm, on the 
other hand, allows to stop and start again the sam- 
pling procedure, in case of charge conservation failure, 
when a number of random extractions smaller than K 
has been performed by checking the conservation of the 
single abelian charges during the sampling. 

The only exception to the above sampling procedure 
is represented by the heavy flavored hadron sampling 
step, required when a heavy flavored cluster must be 
hadronized. For these objects, which are supposed to 
hadronize into channels containing one heavy hadron 
and a set of light particles, the heavy hadron is sampled 
with probability Ph of the form 



P h oc e~ 



(13) 



with a > and where m,h is the hadron mass. This par- 
ticular choice for the heavy flavored hadron sampling 
function is justified by two reasons: on one side it allows 
to avoid the strong undersampling which would follow 
from the usage of the standard sampling method (Eq. 
[9| also for these hadronic species. On the other side, 
the above function allows to approximativcly take into 
account the heavy hadron relative production probabil- 
ities due to their mass differences. The remaining part 
of the channel is then randomly chosen using the same 
algorithm applied for the hadronization of light clus- 
ters. This modification in the sampling function for the 



heavy flavored clusters allows to obtain a more efficient 
channel generation algorithm with respect to the stan- 
dard procedure previously described. 

Given a hadronization channel composition, a kine- 
matical configuration must be generated: in the present 
work, a sampling algorithm inspired by the multi-parti- 
cle decay phase space integration described in 9j^| is 
used. The starting point to obtain the kinematical con- 
figuration and the corresponding weight, for a decay 
channel of a cluster of mass M, containing N particles 
with masses mi, m2, mjv, is the phase space integral 
of Eq. [T| 



PS(M, m 1 ,m 2 , m N ) 



= [f[d 3 Pi s* (p-jrp) 

J »=i V t=i / 



(14) 



As discussed in [5] , the following condition holds for the 
relativistic iV-body phase space element: 



&$n(P\Px,P2, -,Pn) = d$N-j+i(P; Q,Pj+i, ...,pn)x 
d<MQ;pi,p 2 ,...,Pj)(27T) 3 dQ 2 

(15) 

for 2 < J < N, where 



(N \ N 
n 
»=i / i=i 



d 3 pi 



{2^f2E l 
(16) 



and 



The non- relativistic phase space integral of Eq. [14] can 
be rewritten making use of the relativistic phase space 
element of Eq. [16] as follows: 



PS(M,mi,iri2, — ,mjv) 



N ,3 / N \ N 



JV 

/ d$ N (P; Pl ,p 2 ,..., PN )Y[(2ir) 3 2E, 

i—1 



(17) 



1 The detailed description of the adopted phase space sam- 
pling algorithm can be found in | 10| . 
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The integration/sampling procedure adopted here is 
based on the property of Eq. |15[ iteratively applied in 
order to deal with 2-body decays only as described by 



PS(M,m 1 ,m 2} ...,m N ) = 



„ N-3 

/ Y[ d^ 2 {q i -p i+1 ,q l+l ){2TTfdq 2 l+1 : 



(18) 



N 



d# 2 (q N -2]PN-i,PN) Y[(2Tr) 3 2Ei, 



with qo — (M, 0) and where the integration limits of 
the qf virtualities are given by 



and Ei is the i-th particle energy in the cluster rest 
frame corresponding to the momentum p.;. 



It follows from Eq. 20 that the needed ./V-body phase 
space configuration can be obtained by sampling iV-1 
solid angles, corresponding to the particle emission di- 
rection in the 2-body decays, and N-2 virtualities. In 
the present work, the solid angles fip. and the qf vir- 
tualities are randomly generated within 47r and within 
the limits defined in Eq. [19] respectively. The parti- 
cle momenta generated with this procedure, initially 
defined in the 2-body decay rest frames, are subse- 
quently boosted to the laboratory frame. Finally, the 
phase space sampling weight wps is given by 



Jl rain 



N 

E 5 



(19) 



2 max ( Pi V2 „ A 

Qi = [vw-v - m i ) 



By solving the 2-body phase space integrals, Eq. \TE\ 
becomes 

PS(M,m,i,m 2 , ...,m N ) = 



N-2 



2 N- 



n / ^ / d °p 



i 



|Pjv-i| 



f ft(0) + «i(0) 



A' 



(20) 



where |pi| is the i-th particle momentum modulus, de- 
termined by energy-momentum conservation, in the i- 
th 2-body decay rest frame, namely 



1/2 



2q t 



(21) 



for i = l, N — 2, and 



w PS 
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(24) 



Defining the hadronization channel sampling weight 
Whc as 



whc 



K 1 



NA 



(2J 3 + 1) V 

(27T) 3 
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(25) 



for light clusters and as 
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NA 



(2 Jj + 1)V 

(27T) 3 



P h n {Nj} 



(26) 



for heavy flavored clusters, where the hadronization 
channel sampling functions II^ N .j and Ph are defined 
in Eqs.[9|and[T3| respectively, the total weight w corres- 
ponding to a single cluster hadronization event reads 



(qn-2 - ( m N-i + m N ) 2 ^ (q 2 N _ 2 - (m N -i - m N ) 2 ^j 



1/2 



WHC'Wps 

Q 



(27) 



2q 



N-2 



(22) 



for i = N — 1 and i = N. Moreover, the energy compo- 
nents pi(0) and <7i(0) are given by 



ft(0) = \/m 2 + |p 4 | 2 
9<(0) 



(23) 



Itf + IP*I 2 



3 Comparison with the Cluster Model 

Although the SHM is based on the general idea of clus- 
ter hadronization, it has peculiar differences with re- 
spect to other hadronic event generators implementing 
Cluster Model hadronization such as Herwig++ [IT] and 
Sherpa [T^], besides Herwig6510. Let us now discuss 
the differences between these approaches more in detail 
starting from the Herwig case. 
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The most striking difference is concerned with the 
generation of the channel multiplicity: whilst the Clus- 
ter Model forces clusters to decay into hadron pairs, 
in SHM the final multiplicity is an outcome of the clus- 
ter hadronization process. Particularly, the relative pro- 
duction rate of A^-body channels is determined by the 
cluster's finite volume. This happens because the SHM 
uses the proper measure of the multihadronic phase 
space [TlfT3] with a cluster's proper volume V acting 
as a coupling constant for the cluster decay according 
to Eq. [T] Conversely, the Cluster Model uses invariant 
momentum space measure d 3 p/2E without extra fac- 
tors and therefore the final multiplicity is to be fixed 
otherwise (see extensive discussion in |13jV 

Indeed, in the Herwig event generator, a hadroniz- 
ing cluster can be transformed in three different ways 
depending on its mass: 

— splitting into a lighter cluster pair (see discussion in 



Sec. 5.2) 



— decay into a hadron pair; 

— transformation into the lightest hadron with the 
same flavor composition of the cluster. 

Here we point out a further difference, even for the two- 
body decays: given a cluster with flavor composition 
fif2, in the Herwig Cluster Model the decay hadron 
pair will have the flavor composition f]ff and ff 2 re- 
spectively, with the hadrons sharing the components 
of the ff parton pair extracted from the vacuum to per- 
form the cluster decay in this framework. On the other 
hand, in the SHM only abelian charges conservation is 
required, implying the inclusion of a larger set of chan- 
nels. 

With respect to the general aforementioned Herwig 
hadronization scheme, Sherpa implementation of the 
Cluster Model features additional differences with re- 
spect to SHM. Indeed, relevant changes have been intro- 
duced in implementing the Cluster Model in the event 
generator [T3], both for cluster production and decay: 
inclusion of soft color reconnection during cluster for- 
mation and decay; inclusion of the parton/hadron spin 
degree of freedom in the transition probability calcu- 
lation; a refinement of the single hadron cluster trans- 
formation. Particularly, in the Sherpa case all clusters 
with mass falling into the so-called hadronic regime are 
transformed into single hadrons with the same flavor 
compositions of the hadronizing clusters. Each specific 
hadron is then randomly chosen among those with a 
mass smaller than the one of the cluster itself, instead 
of always choosing the lightest one as in the Herwig 
case. 



4 Multi-cluster hadronization and event weight 

Within cluster hadronization models, high energy e + e _ 
collisions (say y/s > 4 GeV) always involve multiple 
cluster production. For instance, the Cluster Model [5 
implemented in Herwig6510 event generator [4], pro- 
vides the formation of color connected parton-antiparton 
pairs at the end of the perturbative QCD shower, gi- 
ving rise to massive colorless clusters. As discussed in 
Sec. [3j those clusters are thereafter decayed into pairs 
of massive hadrons, the probability of the specific decay 
channel being determined by the two-body phase space 
availability and by the spin multiplicities of the involved 
hadrons. Indeed, in the present work, the multi-cluster 
system formed as a result of the hard scattering pro- 
cess, followed by the perturbative showering, is taken 
from the Herwig6510 event generator and used as in- 
put for the SHM. Specifically, the general simulation 
setup is as follows: the Herwig6510 event generator 
is used to simulate the full pre-hadronization process, 
from initial state radiation emission to cluster produc- 
tion. The Herwig6510 preliminary cluster splitting pro- 
cedure described in £Q , which follows the cluster forma- 
tion step, is also included in the event simulation (the 
details of the cluster splitting procedure and its impact 
on the SHM predictions are discussed in Sec. 5.2 1. At 



that stage, the standard hadronization algorithm is re- 
placed by the SHM one, which takes care of hadronizing 
all produced clusters; finally, the formed hadrons are 
returned to the external generator (i.e. Herwig6510) 
which takes care of decaying of the unstable particles. 

The above hadronization procedure is applied to 
each cluster independently, thus, for each simulated col- 
lision, the hadronization event weight W reads: 



N c 

n 



w„ n w 



N c 



HC W PS 



(28) 



where Nq is the number of clusters in the event and Wi 
the i-th cluster hadronization weight defined in Eq. |27| 
Moreover, w l HC and Wp S are the hadronization channel 
and phase space sampling weights of the i-th cluster 



(defined in Eqs. 24 25] and 26 1 and Qi its partition 



function (defined in Eq. [3| . 

This method in fact involves a complication, related 
to the calculation of the hadronization event weight of 
Eq. [28j Specifically, it requires the prior knowledge of 
each cluster's microcanonical partition function. If we 
had only single cluster event with fixed mass and char- 
ges, its microcanonical partition function would be an 
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irrelevant constant factor cancelling out when calculat- 
ing mean values. On the other hand, for a multi-cluster 
environment, the value of the product of microcanoni- 
cal partition functions in Eq. [28] is not constant, in fact 
it is a function of the masses and the charges of the par- 
ticular set of clusters. Therefore, we need to know the 
specific numerical value of the microcanonical partition 
functions of the clusters produced in a single event in 
order to correctly normalize its weight. 

Unfortunately, computing the microcanonical parti- 
tion function of a cluster of given mass and set of char- 
ges is a CPU time demanding task that cannot be af- 
forded at the event generation time. To solve this prob- 
lem, we decided to pre-calculate microcanonical parti- 
tion functions for a discrete set of masses (and charges) 
and store their values in a look-up table to be used 
at event generation time. The number of pre-calculated 
functions has been determined on the basis of the obser- 
vation that, in about 10 6 e + e~ collision events at the Z 
peak, around 150 different cluster charge configurations 
(including baryon number, strangeness, electric charge, 
charm and beauty) occurred. For each of these charge 
configurations, a set of microcanonical partition func- 
tions has been calculated, each set including different 
cluster mass and different free parameters (73 and p) 
values. Specifically, for each charge configuration, a grid 
in the mass-/3-7s space has been defined as follows: 

- p £ [0.20; 0.50] GeV/fm 3 with Ap = 0.05 GeV/fm 3 ; 

- 7s € [0.50; 1.00] with Z\ 7s = 0.05; 

- Cluster mass M £ [277V; 10 GeV] with: 

- AM = 0.1 GeV for M < 3 GeV; 

- AM = 0.5 GeV otherwise; 

and for each grid point an fl calculated. Then, du- 
ring event generation, the partition function values of 
the single clusters, for their specific charge configura- 
tion/mass and for the chosen SHM free parameter val- 
ues, are determined by means of a linear interpolation 
between the nearest points in the grid. 

As it can be seen, a non uniform grid structure along 
the cluster mass direction has been adopted: this choice 
was motivated by the non-smooth dependence of the 
microcanonical partition function on the mass at low 
mass values. Indeed, 3 GeV as mass value separating the 
two differently discretized regions proved to be a good 
trade-off between the need of good accuracy and the 
computational cost in getting all grid points calculated. 

A final remark concerning the calculation of the par- 
tition function set is in order: for a given mass, charge 
configuration and 7s and p parameters values, the cor- 
responding microcanonical partition function value is 



determined by the set of hadrons used in the hadroniza- 
tion process and by their physical properties. It follows 
that the set of partition functions needs to be rebuilt 
only when changes are introduced in the hadron set, 
such as modifications to the hadron physical properties 
or to the hadron list. 

The details of the numerical methods used for the 
pre-calculation of microcanonical partition function grid 
can be found in |10j . It is worth mentioning here that 
we have used the phase space integration optimized al- 
gorithm described in [8TT5"]. 

5 Preliminary tests 

In this section the results obtained with the SHM code, 
interfaced to Herwig6510, for e + e~ collisions at 91.2 
GeV center of mass energy, will be presented and dis- 
cussed. While a rigorous tuning of the SHM free para- 
meters is needed to assess the goodness of the model, 
the results presented here, obtained with an approxi- 
mate parameter adjustment, already show a good ac- 
curacy level. Moreover, in view of a global tuning of 
both SHM and Herwig6510 free parameters, a prelimi- 
nary analysis of the dependency of SHM predictions on 
some Herwig6510 parameters involved in the cluster 
formation step is presented. 

5.1 Simulation setup 

As already discussed, the Herwig6510 event generator 
has been used for the present work as external code for 
the collision simulation, replacing its standard hadroni- 
zation routines with the SHM Monte-Carlo module. As 
input clusters for the microcanonical hadronization, the 
clusters normally produced by Herwig6510 have been 
used, except for two modifications described below. The 
set of hadrons produced during each microcanonical ha- 
dronization process is then processed by Herwig6510, 
which performs the unstable particle decays. 

The modifications introduced in the standard clus- 
ter sets produced by Herwig6510 during the event si- 
mulation, and above mentioned, are the following: 

— Cluster merging: the predictions of the SHM show 
a dependency on the input cluster mass spectrum, 
which in its turn depends on the Herwig6510 setup. 
With the final aim of obtaining a better agreement 
between the SHM predictions and the experimental 
data, it is useful to have the possibility to modify 
the properties of the cluster mass spectrum. In the 
present work, the above modification is introduced 
by means of a cluster merging procedure, which, 



working after the Herwig6510 cluster production 
step, allows to control their mass spectrum with- 
out any changes in the external generator standard 
setup. Starting from a set of Herwig6510 clusters, 
the result of the merging procedure is determined by 
the value of a free parameter (Mq) representing the 
minimum allowed cluster mass: the set of incoming 
clusters is analyzed in order to check the presence 
of objects with mass under the given threshold. If 
one or more light clusters are found, an iterative fu- 
sion process is activated, which merges the cluster 
pairs containing at least one light object into heav- 
ier clusters, repeating this procedure until no light 
clusters remain. When more than one combination 
for a light cluster is possible the merged pair is the 
one with the smallest invariant mass 



m(i,j) 



(Pi +Pj) 



(29) 



where p%(pj) is the 4-momentum of the i-th(j-th) 
cluster. It must be noted that a maximum mass 
value of 10 GeV is allowed for the new clusters pro- 
duced during the merging procedure, a condition re- 
quired to control the broadening of the cluster mass 
spectrum towards large mass values and the conse- 
quent quick increase in the computational cost of 
the corresponding partition function grid calcula- 
tion. Because of the above condition, when a clus- 
ter pair is chosen for the merging procedure, the 
corresponding new cluster is created only if the pair 



invariant mass (Eq. 29) is smaller than 10 GeV: in 



a limited number of cases this condition can pre- 
vent the cluster merging procedure from being fully 
completed, with the result that clusters with mass 
below the chosen Mq value can still be included in 
the hadronization procedure. 

Baryonic clusters: the second change to the stan- 
dard Herwig6510 cluster properties is related to the 
production of baryonic clusters, disabled in a stan- 
dard run of this generator, and is obtained by means 
of a proper setting of its QDIQK parameter value. 
This parameter represents the maximum scale at 
which gluons can be (non-perturbatively) split into 
diquarks during the QCD shower process and is set 
to zero in the default Herwig6510 configuration, 
thus switching off the possibility of diquarks produc- 
tion. For the present work the value 2m c — 2 x 1.55 
GeV has been chosen for the QDIQK parameter, with 
the result of activating the possibility of gluon split- 
ting into light flavored diquark-antidiquark pairs. 
During the subsequent Herwig6510 cluster building 
step, each produced diquark (antidiquark) is cou- 
pled to the corresponding color connected quark 



(antiquark) to form a colorless baryonic cluster. The 
need for baryonic clusters is strictly related to the 
correct baryon production during the microcano- 
nical hadronization process: because of the bary- 
onic charge conservation, baryons can be obtained 
from the standard Herwig6510 non-baryonic clus- 
ters only as baryon-antibaryon pairs. However, these 
configurations are strongly suppressed because of 
phase space availability reasons. The introduction 
of baryonic clusters, on the other side, allows to ob- 
tain the production of baryonic final states in less 
restrictive phase space conditions. 

In Fig. [I] the effects of the discussed cluster mer- 
ging procedure on the cluster invariant mass spectrum 
are shown, in particular for clusters produced in 10 6 
c+c~ — > dd collisions at 91.2 GeV. Primary clusters are 
the standard Herwig6510 clusters, except for the acti- 
vation of baryonic cluster production, and their mass 
spectrum is compared with the one of clusters obtained 
with the merging procedure for Mq = 1.6 GeV: the pre- 
sence of the cluster minimum mass selection is evident, 
as well as of a broadening of the mass distribution. At 
the same time, it can be seen that a residual fraction of 
clusters, about 1%, have a mass below the Mq value: 
the existence of these objects after the merging proce- 
dure is due to the maximum allowed cluster mass value 
of the merging procedure which, as already discussed, in 
a limited number of cases prevents the cluster merging 
from being completed. The mass distribution of the sec- 
ondary clusters, namely the ones provided in output by 
the merging procedure, shows also an inflection point 
around 3.5 GeV. As it can be seen in Fig. [TJ this mass 
value is the maximum allowed for primary clusters, a 
condition responsible for the observed inflection point: 
in fact, only new clusters produced by the merging pro- 
cedure contributes to the mass distribution for values 
larger than the above limit, while for smaller mass val- 
ues contributions from the primary cluster mass distri- 
bution are also present. These differences in the sec- 
ondary cluster mass distribution composition, in the 
two mass ranges separated by the primary cluster max- 
imum mass value, is the origin of the behavior change 
in the mass distribution. 

Finally, we point out that in the hadronization step, 
among the light flavored hadron states, only those with 
mass less or equal to 1.8 GeV were included. This choice 
was motivated by the opportunity of comparing multi- 
plicities with previous results obtained with the SHM 
where this cutoff was used. It should be noted, any- 
how, that the production of light flavored states with 
mass larger than 1.8 GeV is negligible for essentially all 
observables. 
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Fig. 1 Cluster invariant mass distribution for 10 6 
e+e - — > dd events at 91.2 GeV: comparison among primary 
clusters and clusters obtained with the merging procedure 
described in the text. 



5.2 Study of Herwig cluster splitting 

As already stated, the Herwig6510 cluster splitting pro- 
cedure described in 0] is included in the present event 
simulation setup. The splitting procedure, which in Her- 
wig6510 precedes the cluster decay step, operates as 
follows: a cluster of mass M f j , composed of the par- 
tons fi and f2 of mass and toj 2 respectively, is split 
if the following condition holds: 



(M flf2 ) CLPDW > CLMAX CLP ™ + (m fl + m h f 



(30) 



where CLPOW and CLMAX are phenomenological param- 
eters. For each cluster satisfying the above condition a 
flavor pair ff is picked from the vacuum and its com- 
ponents used for the construction of the two daughter 
clusters, whose flavor compositions will be fif and ff2 
respectively. The masses M f j and M f f of the newly 
created clusters are randomly generated according to 
the equations 



M fJ = m fl + (M iA - m fl - m f )r 1 / pspLT 
M fh = mj 2 + (M fl?2 - mj 2 - m f )r 1 / psPLT 



(31) 



where mf is the mass of the generated parton pair com- 
ponents, r is an uniformly distributed random num- 
ber and where PSPLT is a third free parameter of the 
splitting procedure (in fact, PSPLT is a two component 
parameter: one component refers to light flavored and 
charmed clusters and the other to the bottomed ones). 

While a detailed analysis of the interplay between 
SHM free parameters and those belonging to Herwig6510 
initial collision steps is needed in view of a global tun- 
ing, especially concerning the cluster formation and split- 
ting free parameters such as CLPOW, CLMAX and PSPLT, 
a preliminary evaluation of their impact on the SHM 
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Fig. 2 Rapidity with respect to thrust axis normalized distri- 
bution: CLMAX parameter variation effects on SHM predictions 
for e+e - — > dd collisions at 91.2 GeV center of mass energy. 



predictions is here reported. The discussed numerical 
results have been obtained by means of the following 
single variations in the above cluster splitting param- 
eter default values (the default condition PSPLT (2) = 
PSPLT (1) has been used in the present analysis): 



- CLMAX 

- CLPOW 

- PSPLT 



2.35 
1.00 
0.50 



3.35 (default) 
2.00 (default) 
1.00 (default) 



4.35 GeV; 

4.00; 

2.00. 



The analysis has been focused on e + e~ — » dd collisions 
at 91.2 GeV and the following values have been used 
for the SHM hadronization module free parameters: 

- M c = 1.6 GeV; 

- p = 0.45 GeV/fm 3 ; 

- 7s = 0.65. 

The obtained results show that only small effects are in- 
troduced by the considered variations in the CLMAX and 
CLPOW parameter values, as can be seen e.g. for the ra- 
pidity distribution reported in Fig.[2]and|3](for the def- 
inition of the observables, see Appendix). On the other 
side a considerable dependency on the PSPLT parameter 
value is observed in the SHM predictions, as clearly visi- 
ble in Fig.[4]again for event rapidity. The same situation 
is present for single particle observables, as reported for 
example in Figs. [5] - [7] for p° meson scaled energy dis- 
tribution, and for particle multiplicities (Tabs. [I] - [3]). 
Concerning multiplicities, it is worth noting the large 
variations in the total number of charged particles and 
in light meson multiplicities due to the modification of 
the PSPLT parameter, while CLMAX and CLPOW variation 
effects are essentially limited to the baryonic sector. 
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Fig. 3 Rapidity with respect to thrust axis normalized distri- 
bution: CLPOW parameter variation effects on SHM predictions 
for e+e - — > dd collisions at 91.2 GeV center of mass energy. 
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Fig. 6 p scaled energy distribution: CLPOW parameter vari- 
ation effects on SHM predictions for e+e - — > dd collisions at 
91.2 GeV center of mass energy. 
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bution: PSPLT parameter variation effects on SHM predictions 
for e+e - — > dd collisions at 91.2 GeV center of mass energy. 
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Fig. 7 p° scaled energy distribution: PSPLT parameter vari- 
ation effects on SHM predictions for e+e - — > dd collisions at 
91.2 GeV center of mass energy. 
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Fig. 5 p° scaled energy distribution: CLMAX parameter vari- 
ation effects on SHM predictions for e+e - — > dd collisions at 
91.2 GeV center of mass energy. 



5.3 Adjustment of the free parameters 

The free parameters of the SHM code, namely p and 
7s, and the cluster minimum mass parameter Mq in- 
troduced by the merging procedure, need to be rigorou- 
sly tuned to correctly evaluate the performances of the 
code in reproducing the experimental data. For the pre- 
sent work only preliminary tests have been performed, 
using different parameter configurations, in order to ob- 
tain a global overview of the SHM predictivity level in 
a full collision simulation framework. Moreover, an ap- 
proximate adjustment of the SHM code free parameters 
has been realized, through a comparison between the 
SHM predictions and the corresponding experimental 
data: this comparison has been performed considering 
a set of inclusive and exclusive observable distributions 
and mean values. In particular, event shape and single 
particle momentum/energy distributions, single hadron 
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Table 1 Charged particles and hadron multiplicity mean values: CLMAX parameter variation effects on SHM predictions for 
e+e - — > dd collisions at 91.2 GeV center of mass energy. 





Default 


CLMAX = 2.35 GeV 


CLMAX = 4.35 GeV 


Charged 


20.29 ±0.32 


19.93 ± 0.10 


19.64 ± 0.17 


7T° 


10.09 ±0.08 


10.10 ± 0.06 


9.95 ±0.09 


7T + 


8.71 ±0.16 


8.55 ±0.05 


8.46 ±0.08 


V 


1.24 ± 0.03 


1.25 ±0.02 


1.21 ± 0.02 


A++ 


0.081 ± 0.002 


0.0724 ± 0.0007 


0.080 ± 0.002 


2+ 


0.0277 ± 0.0006 


0.0245 ± 0.0003 


0.0229 ± 0.0004 


2- 


0.0243 ± 0.0005 


0.0233 ± 0.0002 


0.0238 ± 0.0004 


2° 


0.0357 ± 0.0008 


0.0266 ± 0.0003 


0.0287 ± 0.0005 



Table 2 Charged particles and hadron multiplicity mean values: CLP0W parameter variation effects on SHM predictions for 
e+e _ — > dd collisions at 91.2 GeV center of mass energy. 





Default 


CLP0W =1.00 


CLP0W = 4.00 


Charged 


20.29 ±0.32 


20.07 ±0.34 


20.00 ±0.12 


7T () 


10.09 ±0.08 


10.05 ±0.14 


10.04 ±0.07 


7T + 


8.71 ±0.16 


8.61 ±0.15 


8.55 ± 0.05 


V 


1.24 ±0.03 


1.28 ± 0.05 


1.22 ± 0.02 


A++ 


0.081 ± 0.002 


0.076 ± 0.002 


0.0864 ± 0.0009 


2+ 


0.0277 ±0.0006 


0.0233 ± 0.0005 


0.0239 ± 0.0003 


2~ 


0.0243 ± 0.0005 


0.0228 ± 0.0005 


0.0246 ± 0.0003 


2° 


0.0357 ± 0.0008 


0.0284 ± 0.0006 


0.0270 ± 0.0003 



Table 3 Charged particles and hadron multiplicity mean values: PSPLT parameter variation effects on SHM predictions for 
e+e _ — > dd collisions at 91.2 GeV center of mass energy. 





Default 


PSPLT = 0.50 


PSPLT = 2.00 


Charged 


20.29 ±0.32 


18.5 ±0.20 


23.1 ± 1.1 


n° 


10.09 ± 0.08 


9.26 ±0.07 


11.6 ±0.5 


7T + 


8.71 ± 0.16 


7.98 ±0.10 


9.9 ± 0.5 


V 


1.24 ± 0.03 


1.14 ±0.03 


1.5± 0.2 


A++ 


0.081 ± 0.002 


0.067 ±0.001 


0.075 ± 0.003 


2+ 


0.0277 ± 0.0006 


0.0205 ± 0.0004 


0.029 ± 0.001 


2~ 


0.0243 ± 0.0005 


0.0203 ± 0.0004 


0.069 ± 0.003 


2° 


0.0357 ± 0.0008 


0.0247 ±0.0005 


0.027 ± 0.001 



mean multiplicities and the charged particle number 
distribution have been included in the analysis. The ref- 
erence experimental data used for the comparison come 
from the measurements of LEP experiments for a center 
of mass energy of 91.2 GeV. The best parameter estima- 
tion has been realized by means of an approximate % 
minimization procedure: a set of explorative runs with 
various free parameter configurations has been consid- 
ered, evaluating for each configuration the discrepancy 
between theoretical predictions and experimental data 
in the form of a global x 2 value. As usual, the param- 
eter best configuration has been identified as the one 
giving the lowest global x 2 value. More in detail, for 
histograms the x 2 value has been computed as the sum 



over channels of the discrepancy between the theoret- 
ical prediction y\ and the corresponding experimental 
data yl, normalized using the theoretical (Monte-Carlo) 
and experimental errors a\ and o\: 

2—1 L 

where i is the histogram channel index and N the num- 
ber of channels. The \ 2 value corresponding to the par- 
ticle multiplicity analysis has been computed in a sim- 
ilar way, summing over the list of considered particles. 
Finally, the global \ 2 value has been computed as the 
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sum over the single distribution/mean value \ 2 contri- 
butions. 

The performed analysis, whose details are reported 
in |10j . shows that the parameter setup corresponding 
to the best global agreement between SHM predictions 
and the corresponding experimental data is: 



M c = 1.6 GeV; 
p = 0.45 GeV/fm 3 ; 
7s = 0.65. 



6 Numerical results 

In this Section, we present the results obtained with the 
SHM code for the above parameter set. The SHM pre- 
dictions are compared to the corresponding predictions 
of Herwig6510 in its release version and to experimen- 
tal data. As already mentioned, the reported results 
refer to e + e~ collisions at the Z peak. It must be noted 
that the chosen observables provide a reliable measure 
for the evaluation of the hadronization model predic- 
tivity, because of their sensitivity to the hadronization 
process. This condition is clear for what concerns exclu- 
sive observables involving a single hadron species, such 
as the tt~^~ moment distribution and mean multiplicity. 
Nevertheless, the same property holds also for the in- 
clusive event shape observables referring to the event 
transverse plane, such as transverse momentum: the 
main contribution to the global hadron distribution in 
the event phase space comes from the hard scattering, 
whose quark emission direction approximatively defines 
the event thrust axis. On the other side, the hadro- 
nic phase space configuration projection on the event 
transverse plane is strongly influenced by the hadroni- 
zation process. Therefore, the analysis of the transverse 
plane related observables allows to effectively evaluate 
the prediction accuracy of the considered hadronization 
model. The results obtained for a subset of these inclu- 
sive and exclusive observables are reported hereafter 
(for the definition of observables, see Appendix). 

— Event shape observables: the comparison between 
the SHM predictions and LEP experimental data 
shows a quite good global agreement for the whole 
set of considered observables (Figs. [8|- [TT]) . However 
the SHM, for the present parameter configuration, 
seems to fail in reproducing the transverse momen- 
tum distribution tail behavior of the experimental 
data (Figs, [s] - [9]) , which is instead correctly pre- 
dicted by Herwig6510. 



Single particle observable^ (Figs. 12 - 27): also in 
this case a good agreement between the SHM pre- 
dictions and LEP data is present for a large set of 
the considered single particle observables. Neverthe- 
less, a wrong behavior in the SHM predictions can 
be seen m the D° and D* scaled energy distributions 



(Figs. 26 - 27 1: in these cases also Herwig6510's pre- 



dictions show a disagreement with the experimental 
data, however the Cluster Model seems to be able to 
reproduce the general behavior of these data better 
than the SHM. A more detailed analysis regarding 
these last results is reported hereafter. 

Charmed hadrons: the performed analysis, involving 
a complete collision simulation, has given the pos- 
sibility to identify some possible limits of the SHM, 
not observed in previous studies. In particular, this 
is the case of the D° and D* scaled energy distribu- 
tions of Figs . [26] and [27| respectively, where the SHM 
fails in reproducing the experimental distribution 
shapes. In these distributions two distinguishable re- 
gions are present: a low energy region (xe ^$ 0.5), 
whose filling events correspond to the D° and D* 
meson production in b-hadron decay, and a high 
energy region approximatively corresponding to the 
scaled energy distribution of the primary D° and 
D* mesons. While the distribution behavior in these 
two regions is qualitatively well reproduced by the 
Cluster Model predictions, the SHM in the adopted 
simulation framework fails in reproducing the high 
energy region. Even though this anomalous result 
needs to be further investigated, it is understood 
that the underestimation of the charmed mesons 
scaled energy, provided by the SHM, is strictly re- 
lated to the large mean number of particles pro- 
duced during the microcanonical hadronization of 
charmed clusters and to the low energy availabi- 
lity condition which follows. A condition not pre- 
sent in the two body cluster decays performed by 
the Herwig6510 Cluster Model. 



Charged particle multiplicity: in this case (Fig. |28[ ) 
a very good agreement between the theoretical pre- 
dictions of the SHM and the experimental data is 
present. It is worth noting that for this distribution, 
in particular for the distribution tails, the predic- 
tions of the SHM show a better agreement with the 
experimental data with respect to Herwig6510's re- 
sults. 

Particle multiplicities: the comparison on the mean 
particle multiplicities (Tabs. |4]-|6|, for this preli- 



2 The theoretical distributions have been normalized to the 
integral value of the experimental ones. 
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minary analysis, shows a global agreement of the 
microcanonical predictions with respect to the expe- 
rimental data within 3a in the 75.5% of the consi- 
dered cases. This percentage becomes 100% if only 
charmed and bottomed hadrons are considered, for 
which a discrepancy within la is present in the 55.6% 
of the cases. The corresponding Herwig6510 predic- 
tions agree with experimental data within 3a in the 
55.1% of cases for the whole set of particles and in 
the 72.2% of cases if only heavy flavored hadrons are 
considered. More in detail, a different behavior be- 
tween the SHM and the Cluster Model predictions 
can be observed for the light flavored baryon mul- 
tiplicities, with the SHM tending to underestimate 
the baryon yields and the Herwig6510 hadroniza- 
tion model showing the opposite behavior. Finally, 
a 5a discrepancy between the SHM prediction of 
charged particle total number and the data can be 
noted. A better agreement for this quantity can be 
obtained with a more refined tuning of the param- 
eters, probably at the cost of a worsening of the 
agreement between data and predictions of other ob- 
servables (e.g. event shape distributions and single 
particle multiplicities) . This discrepancy will be fur- 
ther investigated in the next SHM parameter tun- 
ing. 



at high momentum. While this anomaly needs to be 
further investigated, it is already clear how its origin 
lies in the larger mean number of particles produced 
by statistical microcanonical hadronization of charmed 
clusters as compared with standard Herwig procedure. 

These problems will be further investigated, and a 
global fine tuning of the free parameters performed. 
We envisage an extension of this hadronization Monte- 
Carlo code to high energy pp and pp collisions, with the 
final goal of achieving a global assessment of the Statis- 
tical Hadronization Model in elementary high energy 
collisions. A first public release of this code with the 
needed interfaces for its usage with Herwig6510 is forth- 
coming. 
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7 Conclusions and outlook 

We have presented a Monte-Carlo implementation of 
Statistical Hadronization Model based on the micro- 
canonical hadronization of massive clusters. The de- 
veloped hadronization code has been interfaced to the 
Herwig event generator providing the massive clusters 
to be hadronized. With this method, a full description 
of an e + e _ collision at high energy can be obtained with 
statistical model based hadronization. 

The results obtained in c + c~ collisions at the Z peak 
have been presented and discussed. While a final and 
rigorous tuning of the hadronization module free pa- 
rameters has not been performed in this work, the pre- 
liminary comparison confirms the good agreement of 
the SHM predictions with the data, already observed 
in previous studies. Specifically, the agreement between 
measured and predicted abundances of hadronic species 
is confirmed, particularly in the heavy flavor sector. 
A good agreement, between SHM predictions and ex- 
perimental data, is also found for a set of hadroniza- 
tion relevant event shape variables and single particle 
energy /momentum distributions, with the only excep- 
tion of the charmed meson scaled energy, where a con- 
sistent discrepancy is observed. Indeed, the SHM seems 
to fail in correctly predicting the momentum spectrum 
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Fig. 8 Transverse momentum (in) with respect to thrust axis 
normalized distribution: comparison among SHM predictions, 
Herwig6510 predictions and DELPHI data |16| . 
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Fig. 11 Thrust normalized distribution: comparison among 
SHM predictions, Herwig6510 predictions and DELPHI data 
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Fig. 9 Transverse momentum (out) with respect to thrust 
axis normalized distribution: comparison among SHM predic- 
tions, Herwig6510 predictions and DELPHI data | 16 |. 



10' 



10' 



> 



-§ io- 2 

ib 



10" 



la- 



's. 



O -O LEP data 
A -A Herwig6510 
• SHM 



10" 10 1 
|p| (GeV) 

Fig. 12 7r^ 3-momentum module distribution: comparison 
among SHM predictions, Herwig6510 predictions and OPAL 
data El. 



■ •< A 






1 














* 
















G LEP data 
A -A Hcrwig6510 
• • SHM 








f 



3 



Fig. 10 Rapidity with respect to thrust axis normalized dis- 
tribution: comparison among SHM predictions, Herwig6510 
predictions and DELPHI data | 16 |. 
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tion: comparison among SHM predictions, Herwig6510 pre- 
dictions and OPAL data \TT\ . 
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Fig. 26 D scaled energy distribution: comparison among 
SHM predictions, Herwig6510 predictions and DELPHI data 

EH. 



io- 



10 J 

Jj 10° 

~fcT 

ib 



■10- 



10" 



10" 



»-"••-•■■* 



G -O LEP data 
A -A Hcrwig6510 
• SHM 



0.2 



0.4 0.G 



0.8 



1.0 



Fig. 27 D* scaled energy distribution: comparison among 
SHM predictions, Herwig6510 predictions and ALEPH data 
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Fig. 28 Charged particle number probability distribution: 
comparison among SHM predictions, Herwig6510 predictions 
and OPAL data E3] . 



Appendix: observable definitions 



The thrust T is defined as: 



N 




where n is a unit vector along the thrust axis, N is 
the number of particles and p^ is the i-th particle 
3-momentum. Thrust major M and minor m are si- 
milarly defined, replacing n with iim (perpendicular 
to n) and with n m = njf x n respectively. 
The rapidity with respect to the thrust axis yx is 
defined as: 



Vt 



\ log 



E 



Pt 



E 



Pt 



where E and px are the particle energy and 3-mo- 
mentum projection along the trust axis respectively. 
The in and out components of the transverse mo- 
mentum with respect to thrust axis, p™ and p^" 4 , 
are defined as: 



Pt ~ 

out 

Pt 



P ■ n M 

= P • n m 



where p is the particle 3-momentum and n^/ (n m ) 
is the thrust major (minor) axis previously defined. 
The scaled energy and momentum xe, x p and £ p 
are defined as: 



x E 



2E 
— log 



{.Xp) , 



where E and p are the particle energy and 3-momen- 
tum respectively and y/s is the collision center of 
mass energy. 
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Table 4 Charged particles, photon and hadron multiplicity mean values for e+e - — > qq collisions, with q = u, d,s,c,b: 
comparison among SHM predictions, Herwig6510 predictions and LEP data 9,24. at 91.2 GeV center of mass energy. The last 
two columns contain the discrepancies, measured in standard deviations, for SHM and Herwig6510 predictions with respect to 
experimental data. 





SHM 


Herwig6510 


LEP data 


^SHM 


^Heruig6510 


Charged 


22.34 ±0.27 


20.45 


20.76 ±0.16 


5.00 


-1.94 


7 


23.64 ±0.33 


20.11 


20.97 ± 1.17 


2.20 


-0.74 


7T° 


10.98 ±0.13 


9.56 


9.61 ±0.29 


4.33 


-0.19 


7T+ 


9.33 ±0.15 


8.16 


8.50 ± 0.10 


4.55 


-3.39 


n 


1.18 ±0.04 


0.63 


1.059 ± 0.086 


1.26 


-5.05 


p + 


1.16 ±0.03 


0.97 


1.20 ± 0.22 


-0.18 


-1.06 


p° 


1.42 ±0.04 


1.00 


1.40 ± 0.13 


0.13 


-3.06 


UJ 


1.29 ±0.03 


0.97 


1.024 ± 0.059 


4.10 


-0.87 


n' 


0.13 ±0.01 


0.10 


0.166 ± 0.047 


-0.67 


-1.40 


/o(980) 


0.12 ±0.01 


0.01 


0.1555 ± 0.0085 


-2.43 


-17.2 


a+ 


0.12 ±0.01 


0.01 


0.135 ±0.054 


-0.26 


-2.30 




0.167 ±0.007 


0.1278 


0.0977 ± 0.0058 


7.16 


5.18 


h 


0.17 ±0.01 


0.169 


0.188 ± 0.020 


-0.74 


-0.97 




0.081 ±0.005 


0.072 


0.165 ±0.051 


-1.63 


-1.82 


fi 


0.019 ±0.002 


0.012 


0.0120 ± 0.0058 


1.19 


0.017 


K+ 


1.12 ± 0.02 


1.05 


1.127 ± 0.026 


-0.31 


-2.91 


K° 


1.07 ±0.11 


0.942 


1.0376 ± 0.0096 


0.30 


-9.97 


K* + 


0.34 ±0.05 


0.273 


0.357 ± 0.022 


-0.37 


-3.81 


K*° 


0.33 ±0.02 


0.274 


0.370 ± 0.013 


-1.92 


-7.40 


Kf 


0.031 ± 0.004 


0.0361 


0.036 ± 0.012 


-0.38 


0.006 


P 


0.45 ±0.02 


0.762 


0.519 ± 0.018 


-2.79 


13.5 


A++ 


0.069 ±0.003 


0.148 


0.044 ± 0.017 


1.47 


6.13 


A 


0.128 ±0.004 


0.322 


0.1943 ± 0.0038 


-11.43 


33.6 


E+ 


0.0253 ± 0.0009 


0.0667 


0.0535 ± 0.0052 


-5.35 


2.55 


E~ 


0.0233 ± 0.0008 


0.0548 


0.0410 ± 0.0037 


-4.68 


3.74 


E° 


0.034 ±0.001 


0.0450 


0.0389 ± 0.0041 


-1.23 


2.71 


E*+ 


0.0176 ± 0.0006 


0.0551 


0.0118 ± 0.0011 


4.62 


39.4 




(4.0 ± 0.1) x 10" 3 


0.0391 


0.01319 ± 0.00050 


-17.7 


51.8 


E*° 


(2.12 ± 0.07) x 10" 3 


1.84 x 10" 2 


(2.89 ±0.50) x 10" 3 


-1.53 


30.9 


n 


(1.09 ± 0.04) x 10" 4 


4.94 x 10 -3 


(6.2 ± 1.0) x 10" 4 


-5.11 


43.2 


n 


0.51 ±0.02 


0.683 


0.991 ± 0.054 


-8.49 


-5.70 



Table 5 Mean values of charmed hadron multiplicities for e+e - — » cc collisions: comparison among SHM predictions, 
Herwig6510 predictions and LEP data |24| at 91.2 GeV center of mass energy. The last two columns contain the discrepancies, 
measured in standard deviations, for SHM and Herwig6510 predictions with respect to experimental data. 





SHM 


Herwig6510 


LEP data 


^SHM 


■^Herwig6510 


D+ 


0.22 ± 0.02 


0.287 


0.238 ±0.024 


— 0.53 


2.03 


D° 


0.54 ±0.04 


0.577 


0.559 ±0.022 


— 0.33 


0.82 


D s 


0.110 ±0.009 


0.112 


0.116 ±0.036 


-0.18 


-0.12 


D* + 


0.19 ±0.02 


0.207 


0.2377 ±0.0098 


-2.40 


-3.11 


D*° 


0.23 ±0.02 


0.210 


0.218 ±0.071 


0.11 


-0.11 


D? 


0.015 ±0.001 


0.022 


0.0173 ±0.0039 


-0.46 


1.30 


D*° 


0.033 ±0.003 


0.030 


0.0484 ± 0.008 


-1.85 


-2.27 


D* 


0.072 ± 0.006 


0.036 


0.069 ±0.026 


0.12 


-1.26 


D Bl 


0.0053 ± 0.0004 


0.0044 


0.0106 ±0.0025 


-2.11 


-2.50 




0.0038 ± 0.0003 


0.0059 


0.0140 ± 0.0062 


-1.65 


-1.31 


A c 


0.13 ±0.01 


0.036 


0.079 ±0.022 


2.19 


-1.94 
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Table 6 Mean values of bottomed hadron multiplicities for e+e - — > bb collisions: comparison among SHM predictions, 
Herwig6510 predictions and LEP data |24| at 91.2 GeV center of mass energy. The last two columns contain the discrepancies, 
measured in standard deviations, for SHM and Herwig6510 predictions with respect to experimental data. 





SHM 


Herwig6510 


LEP data 


^SHM 


^Heruig6510 


(B° +B+)/2 


0.411 ±0.005 


0.4471 


0.399 ±0.011 


1.03 


4.37 


B a 


0.105 ±0.001 


0.108 


0.098 ±0.012 


0.61 


0.82 


B* /Buds 


0.69 ±0.01 


0.424 


0.749 ± 0.040 


-1.53 


-8.13 


B** 


0.183 ±0.002 


0.143 


0.180 ±0.025 


0.11 


-1.48 


(B 2 * +Bi) 


0.121 ±0.001 


0.094 


0.09 ± 0.018 


1.73 


0.19 


B* s2 


0.00776 ± 0.00009 


2 x 10" 8 


0.0093 ± 0.0024 


-0.64 


-3.88 


b-baryon 


0.110 ±0.001 





0.103 ±0.018 


0.37 


-5.72 
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